Inverse turbulent cascades and confer mally invariant curves 
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We offer a new example of conformal invariance far from equilibrium — the inverse cascade of 
Surface Quasi-Geostrophic (SQG) turbulence. We show that temperature isolines are statistically 
equivalent to curves that can be mapped into a one-dimensional Brownian walk (called Schramm- 
Loewner Evolution or SLE re ). The diffusivity is close to re = 4, that is iso-temperature curves 
belong to the same universality class as domain walls in the 0(2) spin model. Several statistics 
of temperature clusters and isolines are measured and shown to be consistent with the theoretical 
expectations for such a spin system at criticality. We also show that the direct cascade in two- 
dimensional Navier-Stokes turbulence is not conformal invariant. The emerging picture is that 
conformal invariance may be expected for inverse turbulent cascades of strongly interacting systems. 

PACS numbers: 47.27.-i 



To identify underlying symmetries is a central problem 
of the statistical physics of infinite-dimensional strongly 
fluctuating systems. Turbulence is a state of such a sys- 
tem which is deviated far from equilibrium and is accom- 
panied by dissipation. Excitation and dissipation usu- 
ally break symmetries such as scale invariance, isotropy 
and time reversibility. In fully developed turbulence, the 
scales of excitation and dissipation differ strongly and 
are separated by the so-called inertial interval. The main 
fundamental problem of turbulence is how universal is the 
statistics of fluctuations in the inertial interval 0. Are 
symmetries, broken by excitation and dissipation, even- 
tually restored in that range Q? Cascades can be direct 
or inverse depending on whether the integral of motion 
is transferred towards small or large scales, respectively. 
Symmetries broken by excitation (scale invariance and 
isotropy) are generally not restored in direct cascades 
due to the existence of statistically conserved quantities 
P, 0] . On the contrary, symmetries are expected to be 
restored in the inverse cascade where one looks at scales 
much larger than the pumping scale. This is consistent 
with the observation that inverse cascades are scale in- 
variant 0,0,0- Moreover, it has been shown recently 
that the statistics of zero-vorticity lines in the inverse cas- 
cade of two-dimensional (2d) Navier-Stokes turbulence 
display conformal invariance (i.e. local scale invariance), 
revealing an unexpected connection with percolation |j| . 
Such tantalizing results, while awaiting a theory capable 
of explaining them from first principles, pose new ques- 
tions: Are there other turbulent flows that share this 
property and do they belong to the same universality 
class of percolation? In this Letter we answer them by a 
numerical investigation of SQG turbulence. This system 
is relevant for geophysical applications Q and qualita- 
tively similar to Navier-Stokes 2d turbulence. We show 
that zero-temperature isolines are SLE4 curves at large 
scales (in the inverse cascade). Therefore, the isolines 



belong to the same universality class as the trace of a 
harmonic explorer, certain isolines of a Gaussian (free) 
field [l^, interfaces in the 0(2) model and frontiers of 
Fortuin-Kasteleyn clusters in the 4— state Potts model at 
the critical point (for an introduction to SLE and sta- 
tistical models see Refs. [3, llfjj and references therein). 
This connection allows to obtain analytical predictions 
for some characteristic exponents of cluster and loop 
statistics that compare well with numerical results. 

The SQG model describes a rotating stably stratified 
fluid with a uniform potential vorticity || . The temper- 
ature is advected along a surface bounding a constant 
potential vorticity interior, 

d t T + V-VT = kV 2 T + f , (1) 

and determines the velocity v = z x V?/>, ip(x,t) — 
J dy\x — y\~ 1 T(y,t). Without dissipation and forcing 
(x = f = 0) the equations admit two positive-defined 
quadratic invariants Z = J T 2 dx and E = J ipTdx. In 
the presence of a forcing that injects temperature fluc- 
tuations at a scale If, a double cascade develops, akin 
to the one observed in 2d Navier-Stokes turbulence: the 
"energy" E flows upscales whereas the "enstrophy" Z 
goes downscales. Here we focus on the inverse cascade. 
Requiring the energy flux to be scale independent one 
gets the scaling law S r T = T{x + r) — T(x) ~ r H with 
H = (i.e. logarithmic correlation functions). Indeed, 
numerical simulations show that the temperature field in 
the inverse cascade displays a self-similar statistics with 
a scaling compatible with dimensional expectations (see 
HQ EE El H3 and Fig.QJ. 

We now consider the connected regions of like-sign 
temperature (clusters) and their boundaries (loops) - 
see Fig. To guess the cluster statistics one needs the 
knowledge of the scaling properties of the temperature 
field. Indeed, for a self-similar field with Hurst exponent 
H = the fractal dimension of loops is 3/2 [Jg. If one as- 
sumes that such loop ensemble has a conformal invariant 
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scaling limit, it should belong to the same universality 
class as loops in the 0(2) model in the dense phase. By 
exploiting the Coulomb gas representation of the latter 
system (with g — 1, 0) and general scaling arguments 
[ljj it is possible to derive analytically a set of scaling ex- 
ponents associated to cluster and loop statistics. These 
include the fractal dimensions of clusters and loops, the 
power-law exponents for the number of clusters of given 
mass and the number of loops of given length, radius of 
gyration, or area. In Figure|3such statistics are displayed 
for Surface Quasi-Geostrophic turbulence and shown to 
be consistent with that of 0(2) model. 

These results give a strong indication that loops, i.e. 
zero-temperature isolines, might be conformally invari- 
ant objects. If this is the case they should be statistically 
equivalent to SLE K curves with the diffusivity k = 4, as 
it is conjectured for the 0(2) model. To verify directly 
this hypothesis, we proceeded as follows. First, we iden- 
tify putative SLE traces. After having isolated zero-field 
lines, a walker follows them as they explore the upper- 
half plane. The search ends when the walker touches 
the real axis at a distance from the origin larger than 
l T . A sample trace is shown in Figure 01 (a). This selec- 
tion procedure, self-consistent for k < 4, yields a set of 
curves in the half-plane, which are expected in the scal- 
ing limit to converge to chordal SLE joining two points 
on the real axis. Second, we extract the Loewner driv- 
ing function from the trace. To this aim, let us con- 
sider chordal SLE in the upper-half plane H from to 
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FIG. 1: The second-order structure function of temperature 
differences and, in the inset, probability density functions for 
r = 0.02, 0.04, 0.06 compared to a Gaussian density (solid 
line). Data have been obtained by a direct numerical simu- 
lations of Q by a pseudo-spectral code in a fully periodic, 
square domain of size 1 with 4096 2 lattice points. Gaussian 
white-in-time forcing / has correlation length // « 1/200. 
The system is kept in a statistically stationary state by sup- 
plementing Q with a linear damping term — T/r that models 
bottom friction and extracts energy at very large scales l T oc r, 
(l T = 1/20 ^ 1/10 depending on r). 



Xrx,. We parametrize the curve by the dimensionless t, 
not to be confused with time in The equation for 
gt(z) which maps the half- plane minus the trace up to 
t into H itself, is d t g t = 2 / {(p' (g t )[ip(g t ) - £ t ]} where 
ip(z) = x 00 z/(x 00 — z). The equation for g t can be solved 
for a constant £: Gt $(z) — x 00 {rjx 00 (x 00 — z) + [x^^z — 
V? + ^( Xoo - z) 2 (x x - v) 2 ] 1/2 }/{xl(xoo -z) + [xt(z - 
T]) 2 +At(x 0C ~z) 2 (x 00 —r)) 2 ] 1 ^ 2 }wheTeri = (p In this 

case the trace is the semicircle joining rj and Xoc . For a 
generic £(t) we partition the interval [0, T] into subinter- 
vals [t n ,t n -|_i) with to — 0, t^ + i = T where we approxi- 
mate the driving function by the constant = £(t n ) and 
express g t as a composition G tjv _ tjv _ 1 ^ JV _ 1 o • • • oG tli5o . 
It is now possible to extract the driving function from 
a candidate SLE trace, approximated by a sequence of 
points {zq,Z\, . . . ,Zn+x), where zq = and zjv+i = 
Xqo. The first step is to identify the unique semicir- 
cle passing through the points Xoc and Z\ [see Fig. 0] 
(b) for an illustration]. This yields the values for t/q — 
V _1 (^o) = [Rezi Xoo - (Rezi) 2 - (Imzi) 2 ] /(a^ - Rezi) 
and h = (Imzi) 2 a;^/{4[(Rezi - Xoo ) 2 + (Imzi) 2 ] 2 }. The 
map Gtj^o is then applied to the points resulting in a 
new sequence, by one element shorter: z' k = Gt lt £ (zk+i) 
with k = 1 . . . N. The operation is iterated on the new 
subsequence of points until one obtains the full set of tk 
and that gives a piecewise constant approximation of 
the driving function. 

The result of this procedure is an ensemble of £(i) 
whose statistics converge, for I 2 <nt< l 2 to a Gaussian 
process with the variance (£ 2 (t)) = Kt and k — 4± 0.2, as 
shown in Figure |3J We conclude that, within statistical 
errors, zero-temperature isolines in the inverse cascade 
of SQG turbulence are locally SLE4 curves. This applies 
also to other T contours provided that T T rms . More- 
over, in the limit of very large system size where T(r) 
tends to a self-affine field and T rms diverges, we expect 




FIG. 2: Temperature clusters in the inverse cascade of SQG 
turbulence. These are connected domains with positive tem- 
perature. Negative temperature regions are black. 
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all iso-level loops to be statistically equivalent (another 
application of SLE to non-equilibrium systems has been 
found recently for spin glasses 

Remark that Surface Quasi-Geostrophic and Navier- 
Stokes systems belong to a class that is uniquely spec- 
ified by the transport equation JQ) and by a linear, 
scale- invariant, local in time relationship between the ad- 
vected field and the streamfunction ip(x,t) = J dy\x — 
y\ a ~ 2 T(y,t) SQG dynamics corresponds to a = 1, 
Navier-Stokes equation to a — > 2, T being vorticity. The 
larg e-scale limit of the Charney-Hasegawa-Mima model 
[2l| , relevant for atmospheric and plasma turbulence, cor- 
responds to a = —2. Dimensional arguments for the in- 
verse cascade in this class of models give S r T oc r H with 
H = (2 — 2a)/3. This exponent can be used to infer the 
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FIG. 3: Cluster and loop statistics for SQG turbulence, (a) 
The average area M versus the radius of gyration R. (b) The 
length of a loop (blue symbols) and its externally accessi- 
ble perimeter which is obtained by subtraction of fjords with 
necks smaller than If (green symbols) versus 7?. (c) Number 
of clusters of area between M and AM. (d) Number of loops 
of length between L and XL. (e) Number of loops of radius 
between R and XR. (f) Number of loops of area between 
A and XA. In all figures A ~ 1.1. The solid lines are the 
theoretical expectations for the 0(2) model. 



dimension of the contour loops (3— H)/2 = (7+2a)/6 
(for < H < 1 and — 1/2 < a < 1) and thus conjecture 
that they converge to SLE K curves with k = 4(1 + 2a) /3. 

Let us stress that the temperature field in our model 
has non-Gaussian statistics, see and Fig.^ A similar 
non-Gaussian form with logarithmic moments holds for 
the height function built on independently oriented loops 
from the 0(n) model 0,H(j — yet it requires n ^ 2 and 
k ^ 4 (for k = 4 the statistics is Gaussian). Should 
our field belong to this class, the difference n — 4 would 




FIG. 4: (a) A putative SLE trace, (b) A cartoon of the 
algorithm used to extract the driving function from the trace 
(see text). The procedure has been checked on an ensemble of 
self-avoiding loops where it yielded the correct value k = 8/3 
with an uncertainty of less than 5%. 
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FIG. 5: Statistics of the driving function. In the main frame 
is shown the diffusive behavior of £(t). Lower-right inset: 
the diffusion coefficient is k — 4 ± 0.2. Upper-left inset: the 
probability density of £t, rescaled by its variance is Gaussian 
(different symbols correspond to t = 0.02, 0.04, 0.08). 
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be much larger than our 5% margin of error (as can be 
inferred comparing Fig. |U with the results of 20]). It 
remains to be understood how such a non-Gaussian field 
can have isolines with the same statistics as the isolines 
of the Gaussian free field. 

Let us briefly compare our findings with other turbu- 
lent systems. In the direct cascade of 2d Navier-Stokcs 
turbulence, the vorticity field has logarithmic correla- 
tions (H = 0) and is characterized by very weak, if any, 
deviations from self-similarity; vorticity isolines have di- 
mension D = 3/2, just as in the inverse cascade of SQG 
turbulence. However, the similarities end here. The loops 
in the direct cascade, shown in Figure El are not SLE 
curves since they are not even scale invariant as seen 
from the multifractal spectrum in the inset in Fig.HJb). 
Therefore, it appears that inverse cascades are more akin 
to statistical mechanics systems and more appropriate for 
conformal invariance. Another necessary condition may 
be strong nonlinearity since the turbulence of weakly in- 
teracting waves is generally not conformal invariant (ex- 
cept when it has logarithmic correlations in 2d). An- 
other relevant example is a passive scalar in a spatially 
smooth random flow (Batchelor regime), which also has 
logarithmic correlation functions 22]. In this case, cas- 
cade direction depends on the compressibility of the flow 
0- By a straightforward applications of the formulas 
from [2i|, 0] one can show that the four-point correla- 
tions in the Kraichnan model are not conformal invariant 
in either direct or inverse cascade. 




FIG. 6: (a) A typical loop in SQG inverse cas- 
cade. Inset: generalized fractal dimensions D q — (q — 
1) _1 linir^oo In Z q (r)/ In r where Z q (r) = Yl t J -i( r ) q an d the 
sum runs over a set of N(r) non-overlapping r-boxes cover- 
ing the curve, and fJ.%(r) is the mass of the curve inside the 
i-th box. (b) A typical vorticity loop in the direct cascade of 
2d Navier-Stokes turbulence, obtained from a pseudo-spectral 
simulation with 1024 2 lattice points. The spectrum of gener- 
alized dimensions shows a clear dependence on q (multifrac- 
tality). It decreases from the fractal dimension Do w 3/2 to 
reach values D q « 1 for large q, as a result of the existence of 
long, almost one-dimensional segments of the loops. 



To conclude, we have found the second example of con- 
formal invariance in turbulence thus showing that nei- 
ther Navier-Stokes turbulence nor percolation are unique 
in this respect. Both cases correspond to inverse cas- 
cades. In the direct /inverse cascade we study statistics 
at the scales which are respectively smaller /larger than 
the pumping correlation scale. It is thus not surpris- 
ing that the direct cascade is sensitive to the statistics 
of the pumping 0, even when there is scale invari- 
ance, conformal invariance is absent as shown here. In 
inverse cascades, short correlated random force imposes 
some degree of locality and yet conformal invariance is a 
remarkable example of emerging symmetry since our sys- 
tems are dynamically nonlocal and far from equilibrium. 
It remains an open question, here as in the Navier-Stokes 
case, whether conformal invariance extends to some field 
correlation functions, and how to identify candidate pri- 
mary fields upon which a conformal field theory for the 
inverse cascade can be built. 
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